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Abstract 

Background: Mangroves are ecologically important and highly threatened forest communities. Observational and 
genetic evidence has confirmed the long distance dispersal capacity of water-dispersed mangrove seeds, but less is 
known about the relative importance of pollen vs. seed gene flow in connecting populations. We analyzed 980 
Avicennia germinans for 1 1 microsatellite loci and 940 Rhizophora mangle for six microsatellite loci and subsampled 
two non-coding cpDNA regions in order to understand population structure, and gene flow within and among 
four major estuaries on the Caribbean and Pacific coasts of Panama. 

Results: Both species showed similar rates of outcrossing (t= 0.7 in A germinans and 0.8 in R. mangle) and strong 
patterns of spatial genetic structure within estuaries, although A germinans had greater genetic structure in nuclear 
and cpDNA markers (7 demes > 4 demes and Sp= 0.02 > 0.002), and much greater cpDNA diversity {H d = 0.8 > 0.2) 
than R. mangle. The Central American Isthmus serves as an exceptionally strong barrier to gene flow, with high 
levels nuclear {F ST = 0.3-0.5) and plastid {F ST = 0.5-0.8) genetic differentiation observed within each species between 
coasts and no shared cpDNA haplotypes between species on each coast. Finally, evidence of low ratios of pollen to 
seed dispersal (r = -0.6 in A germinans and 7.7 in R. mangle), coupled with the strong observed structure in nuclear 
and plastid DNA among most estuaries, suggests low levels of gene flow in these mangrove species. 

Conclusions: We conclude that gene dispersal in mangroves is usually limited within estuaries and that coastal 
geomorphology and rare long distance dispersal events could also influence levels of structure. 
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Background 

Mangrove communities are composed of phylogenetic- 
ally unrelated species, each adapted in different ways to 
the coastal environment where high salinity, diurnal pat- 
terns of submergence, wave action, and frequent disturb- 
ance create high stress environments [1,2]. Mangrove 
forests have an extended tropical and subtropical geo- 
graphic distribution during the last 40 My [3,4]. They 
are also among the most biologically productive forests 
and provide key ecosystem services such as breeding 
grounds for fish, shrimp, and birds and function to pro- 
tect coasts from tidal surges during hurricanes or 
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extreme weather [5,6]. However, mangrove forests are 
being destroyed and fragmented at alarmingly high rates 
in the tropics due to human development [7-10]. More- 
over, mangrove forests may be particularly susceptible 
to increasing sea levels [11]. The ability of mangroves 
to adapt and regenerate in the face of disturbance 
and changing environmental conditions is therefore 
dependent upon rates of local and long distance disper- 
sal, colonization and genetic diversity found within and 
among populations. 

In most trees, it is generally assumed that rates of 
pollen gene dispersal are greater, and often much 
greater, than rates of gene flow via seed. However, sev- 
eral lines of evidence suggest that mangrove species are 
capable of frequent long distance seed dispersal (LDD) 
on scales of many km that may exceed rates of gene flow 
via pollen. For example, exceptionally high levels of gene 
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flow have been reported between West Africa and South 
America at scales > 6,000 km [12,13]. In addition, mark 
and recapture experiments in ocean currents in the 
Caribbean and in the Gulf of Mexico have revealed high 
dispersal potential for mangrove seeds [14-17]. Direct 
measurements at very local scales within mature man- 
grove forests have shown very short seed dispersal dis- 
tances i.e. less than 10 m in two weeks, [18]. Therefore, 
the extent of gene flow via seed observed among popula- 
tions may be due to the frequency at which seeds en- 
counter large open ocean currents when dispersed from 
more isolated estuaries where water flow is largely due 
to diurnal tidal action. 

A comparison between the extent of pollen vs. seed 
gene flow in mangroves would improve our understand- 
ing of mangrove population structure and reveal the 
contribution of each to observed spatial patterns of gen- 
etic structure [19-21]. In this study, we compare patterns 
of nuclear and plastid genetic structure in the two dom- 
inant mangrove species from the new world, the black 
mangrove {Avicennia germinans, Avicenniaceae) and the 
red mangrove {Rhizophora mangle, Rhizophoraceae), 
across the four main estuaries of Panama within a scale 
up to 300 km along same coastline (Figure 1). The focal 
species currently occur in sympatry on both sides of the 
Central American Isthmus (CAI) but have different histor- 
ical biogeographies. Rhizophora mangle has been present 
in the Neotropics for 40 My while A. germinans has been 
present since 16 My [3,4,22]. The two species have con- 
trasting patterns of life-history traits and current demog- 
raphy. Avicennia germinans has an entomophilous 
pollination system with polyphile (i.e. bees, wasps and flies 
identified as effective pollinators), whereas R. mangle is 
characterized by a simultaneous wind (anemophily) and 
entomophilous pollination, termed ambophilous pollin- 
ation [2,23-26]. Although the two mangrove species have 
perfect flowers, the breeding mechanism seems to be 
different [2]. The protandry observed in A. germinans sug- 
gest a mostly outcrossing breeding system, however self- 
compatibility has not been tested yet in this species [24]. 
In comparison, R. mangle shows a mixed-mating breeding 
system where self-pollination could be more frequent 
[23,27]. Avicennia germinans has small light ovoid crypto- 
viviparous propagules with longevity of four months while 
R. mangle seeds are large, elongate and viviparous with 
longevity of one year [2,28,29]. In addition, while R. 
mangle is the only species from the genus Rhizophora to 
inhabit the Caribbean coast of Panama, on the Pacific side 
it is in sympatry with R. racemosa, generating introgres- 
sive hybrid zones where morphological and genetic dis- 
tinctions between species and hybrids are intricate [13]. 
Finally, populations of A, germinans have a patchy, low- 
density distribution (e.g. 6 stems/km 2 in Bocas del Toro, 
Panama) compared to R. mangle populations, which have 



a more uniform distribution with extremely high densities 
(e.g. 1,544 stems/km 2 in Bocas del Toro, Panama) [30]. 
Previous analyses of A. germinans have reported high 
levels of genetic diversity of this species on the Pacific side 
of CAI, including Panama [12,31,32]. 

Based on these life-history characteristics and previous 
reports of seed dispersal capacity, we expected that: i) 
both species would display low genetic structure or lit- 
tle evidence for isolation by distance (IBD) within estu- 
aries and ii) both species would have a higher rate of 
seed gene flow than pollen gene flow because of their 
capacity for LDD via sea-drift seeds [12,13]. In addition, at 
the species level comparison we expected that Avicennia 
germinans would have greater degree of population 
differentiation than R. mangle because a combination 
of lower seed longevity and viability in ocean water, 
insect pollen dispersal mechanism, and low-density 
distribution [19,33]. 

Results 

Genetic differentiation and mating systems 

Despite differences in many life history characteristics, 
the two species did not show significant differences in 
levels of genetic structure (F ST = 0.32 ± 0.04 SE for A. 
germinans and F S t = 0.40 ± 0.05 SE for R. mangle), 
inbreeding (F IS = 0.20 ± 0.03 SE for A. germinans and 
F IS = 0.13 ± 0.13 SE for R. mangle) or outcrossing rates 
(t = 0.67 ± 0.06 SE for A. germinans and t = 0.77 ± 0.19 
SE for R. mangle) (Kruskal-Wallis ANOVA, P >0.05, 
N=4) (Table 1). Both species showed deviations from 
HWE, but only R. mangle showed evidence of linkage 
disequilibrium. These deviations could be explained by 
deme structure within estuaries that we explained below. 

General patterns of diversity among four estuaries 

Both species showed striking differences in genetic di- 
versity. Avicennia germinans showed considerable vari- 
ation in microsatellite alleles among four estuaries, with 
the highest H e observed in the Caribbean Costa Arriba 
estuary (0.730) and lowest in the Pacific San Miguel Gulf 
(0.459). In contrast, R. mangle showed lower overall 
gene diversity in Caribbean populations (0.349 and 
0.305 for Bocas del Toro and Costa Arriba, respect- 
ively) relative to Pacific estuaries (0.654 and 0.610 for 
Montijo Gulf and San Miguel gulf, respectively), how- 
ever, no significant differences in outcrossing rates 
within species were observed among estuaries (Table 1 
and Additional file 1). The higher gene diversity and 
number of alleles in Pacific estuaries of R. mangle is com- 
plicated by the complex hybridization evident in both 
Pacific estuaries (Figure 1), which serves as a novel source 
of genetic diversity in Pacific populations [13]. 

Chloroplast DNA also revealed striking patterns in di- 
versity and structure in these two species when comparing 
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Figure 1 Population structure for two mangrove species within each estuary when analyzed separately using GENELAND 2.0.12. The 

estuaries under study were Bocas del Toro (a) and Costa Arriba (b) in the Caribbean, and Montijo Gulf (c) and San Miguel Gulf (d) in the Pacific. 
Points represent individuals and their colors represent the assignment of each one of the inferred clusters (K). Admixed individuals that were 
simultaneously assigned to two clusters within each estuary with a probability > 0.5 are in black. Due to the scale of estuary maps, all individuals 
from same transect (i.e. -30 individuals) are overlapping. However, they were usually assigned to the same cluster (i.e. same color in the map) 
with very few exceptions. 



among estuaries. Rhizophora mangle has a single fixed 
haplotype in the Caribbean coast that include Bocas del 
Toro and Costa Arriba estuaries (H d = 0.00), and only two 
haplotypes in the Pacific (H d = 0.44 for Montijo Gulf and 
H d = 0.37 for San Miguel Gulf) for a total of three 



haplotypes. The opposite trend was observed in A. germi- 
nans, which showed remarkably high genetic diversity of 
cpDNA haplotypes with 22 haplotypes observed from 
58 individuals. For this species, the haplotype diversity 
was similar among four estuaries (H d = 0.84 for Bocas 
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Table 1 Nuclear microsatellite F ST and Chloroplast (cpDNA) F ST estimated as G ST in two mangrove species 



Mangrove species Microsatellite 



CpDNA 



F S t Fis Outcrossing level Hamilton & Miller's method 

(SE) (SE) (t) (SE) 



Expected 95% CI 
cpDNA F ST 



Observed 95% CI 
cpDNA G ST 



Ennos 's method 

r = Pollen flow/Seed 
flow 



Avicennio germinons 


0.3181 


0.1950 


0.6736 


0.5550 0.4758-0.6342 


0.3469 


0.1407-0.5531 


-0.64 




(0.0396) 


(0.0335) 


(0.0597) 










Rhizophora mangle 


0.4001 


0.1327 


0.7657 


0.6485 0.5461-0.7509 


0.8504 


0.5888-1.1120 


7.65 




(0.0512) 


(0.1307) 


(0.1885) 











The estimates of G ST were calculed using the Hamilton and Miller's method [21, eq. 10] and the Ennos's method [43, eq. 5a] in Avicennia germinans (black 
mangrove) and Rhizophora mangle (red mangrove) across four Panamanian estuaries. The comparison of F ST , F IS and outcrossing level (t) between A. germinans 
and R. mangle did not show significant differences between species (Kruskall-Wallis ANOVA, P>0.05, N=4). 



del Toro, H d = 0.89 for Costa Arriba, 0.70 for Montijo 
Gulf and 0.73 for San Miguel Gulf). Finally, no cpDNA 
haplotypes were shared between coasts for either spe- 
cies (Figure 2, Additional file 1, Additional file 2 and 
Additonal File 3). 



Hierarchical structure of genetic variation based on 
microsatellites 

We found evidence of similar levels of genetic structure be- 
tween species overall (F S t = 0.32 ± 0.04 SE for A. germinans 
F ST = 0.40 ± 0.05 SE for R mangle, Table 1). In addition, 



Caribbean - 



Bocas del Toro (n=1 5) Costa Arriba (n=1 5) 



Pacific 



Caribbean 



{ 



Pacific - 




Montijo Gulf (n=14) 



San Miguel Gulf 
(n=14) 



Bocas del Toro (n=14) Costa Arriba (n=14) 




Montijo Gulf (n=14) 



San Miguel Gulf 
(n=18) 



Figure 2 Median joining networks and geographic distribution of cpDNA haplotypes across four estuaries in Panama, (a) Avicennia 
germinans (Black mangrove) and (b) Rhizophora mangle (Red mangrove). No haplotypes were shared across the Central American Isthmus. 
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both species showed that majority of microsatellite variation 
is partitioned within estuaries (70% for A. germinans and 
64% for R. mangle) with differences among estuaries 
accounting for most of the remaining variation (30% and 
36%, Table 2). Thus, pairwise comparison of estuaries sepa- 
rated by the ocean showed exceptionally high levels of differ- 
entiation with pairwise F ST ranged between 0.31 to 0.46 in 
A. germinans and between 0.47 to 0.51 in R. mangle. Cor- 
rected estimates of F ST by null allele presence were highly 
similar to non-corrected F S t values putatively harboring null 
alleles (Table 3, Additional file 4). 

In spite of these similar strong structures across 
oceans and estuaries, both species showed differences in 
how genetic diversity is subdivided between estuaries 
from the same coast. Pairwise F ST indicates that popula- 
tion differentiation is lower between estuaries in the 
Caribbean coast than between estuaries in the Pacific 
coast (i.e. F ST = 0.07 and F ST = 0.008 between two Carib- 
bean estuaries and F S t = 0.134 and F S t = 0.08 between 
two Pacific estuaries for A. germinans and R. mangle 
respectively). 

The Bayesian structure analysis supports differences in 
how structure is organized at two levels (2) Between es- 
tuaries from the same coastal line and (3) within estuar- 
ies. Avicennia germinans, with the exception of the 
Costa Arriba estuary, shows evidence of two demes 
within each one of the other three estuaries (Figure 1 
and Additional file 5). In comparison, R. mangle showed 
evidence of three demes within each one of the two 
Pacific estuaries but only one extended population 
across Caribbean coastline without any differentiation 
between Bocas del Toro and Costa Arriba estuaries and 
without any substructure within each one of these two 
estuaries (Figure 1 and Additional file 6). 

At level (3) within estuaries, the two mangrove species 
showed statistically significant negative slopes (b id ) in 
the kinship-distance curves (Figure 3 and Table 4). The 
strength of spatial genetic structure (SGS) measured by 
Sp statistics for the whole estuary is identical if we com- 
pare it with Sp of each deme inferred within the estuary. 



The exception was the Montijo Gulf estuary for both A. 
germinans and R. mangle. In both cases, Sp calculated 
for whole estuary was higher than when Sp was calcu- 
lated for each deme within that estuary. In spite of that, 
we did not find differences in Sp values across seven A. 
germinans s demes distributed in the four estuaries 
(Kruskall-Wallis ANOVA, P = 0.479). However, one of 
the R. mangles demes localized in San Miguel Gulf 
(Sp=0.04 ± SE 0.002) showed a stronger SGS compared 
to only one R. mangle deme shared between estuaries 
from Caribbean coastal line (Sp =0.01 ± SE 0.004) 
(Kruskall-Wallis ANOVA, P = 0.019) (Figure 1, Figure 3 
and Table 4). 

The extent of SGS of two mangrove species was also 
different. Avicennia germinans showed a range of kin- 
ship values almost ten times higher than R. mangle. 
SGS was greater in A. germinans (Sp = 0.0186 ± SE 
0.0026, N demes = 7) than for R. mangle (Sp = 0.0019 ± 
SE 0.0031, N d emes = 4) across demes with N > 50 on scales 
up to 10 km (Kruskall-Wallis ANOVA, P = 0.0000) (Figure 3 
and Table 4). 

Hierarchical structure of genetic variation in cpDNA 

Avicennia germinans has lower genetic structure at the 
cpDNA level (G ST = 0.35) relative to R. mangle (G ST = 0.85) 
(Table 1). However, hierarchical organization of genetic 
structure of the two species at cpDNA shows more similar 
patterns than those inferred by nuclear microsatellites. In A. 
germinans, both Caribbean and Pacific coast samples 
showed strong structure between estuaries, (i.e. F ST = 0.149 
between Bocas del Toro and Costa Arriba estuaries in the 
Caribbean, and F ST = 0.528 between Montijo Gulf and San 
Miguel Gulf in the Pacific). However, some haplotypes were 
shared between estuaries along same coastline but separated 
by -300 km. In contrast, all R. mangle individuals shared a 
single haplotype on the Caribbean coast, while on the Pacific 
side although both estuaries shared the same two haplo- 
types; strong genetic population structure existed between 
Montijo Gulf and San Miguel Gulf (i.e. F ST = 0.353) (Table 2 
and Figure 2). Most of the genetic variation in cpDNA 



Table 2 Analysis of molecular variance (AMOVA) for two mangrove species after 10,000 permutations using ARLEQUIN 
3.5 



Mangrove species 


Microsatellite 






cpDNA 








Variance 


% total 


Fst 


P 


Variance 


% total 


Fst 


P 


Avicennia germinans 


















Among estuaries 


1.113 


30.39 


0.304 


0.000 


2.345 


66.80 


0.668 


0.000 


Within Estuaries 


2.550 


69.61 






1.165 


33.20 






Rhizophora mangle 


















Among estuaries 


0.819 


35.99 


0.399 


0.000 


2.875 


56.23 


0.562 


0.000 


Within Estuaries 


1.229 


64.01 






2.238 


43.77 







Analysis of species Avicennia germinans (black mangrove) and Rhizophora mangle (red mangrove) was performed comparing four estuaries of Panama based on 
microsatellite and cpDNA polymorphism. 
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Table 3 Pairwise genetic structure for two mangrove species across four estuaries in Panama using F ST 







Caribbean 








Pacific 






Avicennia germinans 




Bocas del Toro 




Costa Arriba 




Montijo Gulf 




San Miguel Gulf 


Caribbean 


Bocas del Toro 


- 




0. 149* 




0.817* 




0.721* 




Costa Arriba 


0.068* 

0.063 (0.035 - 


0.099) 


- 




0.740* 




0.658* 


Pacific 


Montijo Gulf 


0.403* 

0.372 (0.298 - 


0.453) 


0.311* 

0.287 (0.214 - 


0.367) 






0.528* 




San Miguel Gulf 


0.461* 




0.369* 




0.134* 










0.442 (0.340 - 


0.549) 


0.361 (0.268 - 


0.460) 


0.140 (0.083 - 


0.201) 




Rhizophora mangle 


















Caribbean 


Bocas del Toro 






0.000 




0.706* 




0.504* 




Costa Arriba 


-0.008 
0.01 1 (0.003 - 


0.019) 






0.706* 




0.504* 


Pacific 


Montijo Gulf 


0.471* 

0.453 (0.327 - 


0.546) 


0.490* 

0.460 (0.334 - 


0.553) 






0.353* 




San Miguel Gulf 


0.494* 




0.512* 




0.078* 










0.470 (0.383 - 


0.555) 


0.479 (0.388 - 


0.564) 


0.071 (0.046 - 


0.104) 





The AMOVA was performed using 10,000 permutations in ARLEQUIN 3.5. It is indicated the F ST values derived from both microsatellite (below the diagonal) and 
cpDNA data (above the diagonal and in italics) for each pairwise comparison across estuaries. The asterkisk indicate P < 0.05, after Bonferroni corrections. 
Underlined values represent the unbiased F ST estimates following the ENA method that correct by the presence of null alleles on F ST estimation (95% CI after 
10,000 replicates). 



genome is partitioned among ocean basins, with 67% for A. 
germinans and 33% of the variation partitioned within estu- 
aries (Table 2 and Additional file 2). Rhizophora mangle 
shows a similar pattern with 56% of the variation partitioned 
among oceans and 44% within estuaries (Table 2). Contrary 
to what was observed across microsatellites, the SGS ana- 
lysis based on cpDNA did not show significant slopes within 
each one of the four Panamanian estuaries analyzed, perhaps 
due to limited sample size (Additional file 3). 

Comparison of pollen and seed migration rates 

In A. germinans the pollen-to-seed gene flow ratio was 
negative, or basically zero (r = -0.64), which indicates 
one to two times higher gene flow via hydrochorious 
seed dispersal than entomophilous dispersed pollen gene 
flow. In comparison, the ratio obtained in R. mangle indi- 
cates approximately seven times higher ambophilous pollen 
gene flow than hydrochorious seed gene flow (r = 7.65). 
However, because the confidence limits of each species 
overlap with the expected value of equal seed and 
pollen flow, we could not reject the null hypothesis of 
wipoiien = J^seed gene flow in either species (Table 1). 

Discussion 

Mangrove communities are critically important ecosys- 
tems that are high in aquatic and terrestrial biological di- 
versity in tropical and subtropical ecosystems worldwide 
[2,5]. Today they are threatened by high rates of anthropo- 
genic disturbance, including habitat destruction, pollution, 



fragmentation, and changes in oceanic and estuarine 
environments due to climate change [9,10,34]. The goal of 
this study was to document spatial genetic structure of 
two dominant Neotropical mangrove species at three 
spatial levels (1) among four estuaries in Panama (2) Be- 
tween two estuaries from same coastal line and (3) within 
each one of the four estuaries. These data provide critical 
information to understanding how genetic diversity is 
structured and maintained within mangrove species and 
communities. 

Both mangrove species showed a strong genetic break 
across the CAI. However, the patterns of diversity 
observed in this study were the opposite of what we had 
expected for both species. We found significant differ- 
ences between estuaries from same coastline and also 
IBD within estuaries in both mangrove species. In 
addition, both mangrove species showed comparable 
outcrossing rates, contrary to other reports on their mat- 
ing systems. Further, A. germinans showed much higher 
levels of genetic diversity, especially in plastid genomes, 
than R. mangle, in spite of the fact that R. mangle 
populations have a more continuous distribution. Finally, 
although there is documented evidence for extreme LDD 
in both mangrove species, our evidence mostly indicates 
restricted gene dispersal overall and largely equivalent 
rates of seed and pollen dispersal. Below, we interpret the 
observed population genetic differences of these species as 
well as the combined effects of species' life histories, gene 
dispersal limitation, and biogeographic history. 
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Figure 3 Spatial autocorrelation of average Kinship coefficients {F u ) against the natural logarithm of spatial distance inside four 
estuaries in Panama. The Kinship-curve for the whole estuary is represented in black. Where GENELAND detected internal substructure, Kinship- 
curves for each genetic pool were calculated when N > 50 (white and grey curves). Dashed lines represent a 95% confidence interval around the 
hypothesis of no genetic structure for the whole estuary based on 10,000 permutations. We generated uneven lags with constant number of 
individuals inside distance classes (N > 100) with > 90% of pairwise relationships among nearest neighbors included within the first interval. 



The ecological imprint on genetic structure: mating 
system and Pollen vs. seed movement 

The two mangrove species showed strong genetic struc- 
ture across four estuaries analyzed, including evidence of 
substructure and IBD within estuaries. Nevertheless, the 
patterns of genetic structure were very different between 
species. Microsatellites revealed lower gene diversity and 
lower genetic structure in R. mangle than A. germinans. 
This contradicts the predictions based upon an outcross- 
ing mating system where assumed outcrossing species, 
such as A. germinans, are expected to have lower genetic 
structure than the mixed-mating R. mangle [33,35]. 



The protandry reported in A. germinans is the main 
evidence that supports an outcrossing mating system in 
this species because the flower- developing mechanism 
makes autogamy unlikely (i.e. within-flower pollination) 
[24]. However, absence of self-pollination (autogamy) is 
not equivalent to self-incompatibility because pollination 
could ocurr between different flowers on the same plant 
(geitonogamy). Other Avicennia species from Indo-West 
Pacific region including A. marina and A. officinalis 
show this pattern [36,37]. There is no direct evidence of 
self-incompatibility in A. germinans that we know of. It 
is possible that the patchy spatial distribution of A. 
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Table 4 Spatial Genetic Structure (SGS) parameters for two mangrove species across four estuaries in Panama 


Mangrove 
species/ Ocean 


Estuary 


N 


bid 

(R 2 ld) 


F A 
(SE) 


Sp 
(SE) 


Avicennia germinans 












Caribbean 


Bocas del Toro 


150 


-0.0217*** (0.0639) 
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0.025 a (0.007) 




Bocas del Toro (P) 


76 


-0.0330*** (0.1038) 


0.108*** (0.023) 


0.037 a (0.009) 




Bocas del Toro (MG ) 


72 


-0.0145*** (0.4050) 


0.104*** (0.016) 


0.01 6 a (0.006) 




Costa Arriba (0) 


241 


0.0081*** (0.0222) 


0.049*** (0.006) 
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Rhizophoro mangle 












Caribbean 


Bocas del Toro and 
Costa Arriba (0) 


All 


— u.uuyj (u.uuoyj 


U.U4o (U.UzzJ 


U.U I U (U.UU4J 


Pacific 


Montijo Gulf 
Montijo Gulf (TG) 


281 

37 


-0.0131*** (0.0241) 


0.054*** (0.011) 


0.0 1 4 a (0.003) 




Montijo Gulf (P) 


233 


-0.0033** (0.0026) 


0.021** (0.006) 


0.003 b (0.001) 




Montijo Gulf (Y) 


10 










San Miguel Gulf 


237 


-0.0276*** (0.0579) 


0.154*** (0.025) 


0.033 a (0.008) 




San Miguel Gulf (R) 


41 










San Miguel Gulf (TG) 


71 


-0.0000108 


0.066*** (0.019) 


0.004 a (0.002) 




San Miguel Gulf (LL) 


125 


-0.00002175 


0.036** (0.027) 


0.008 a (0.011) 



It is indicated the slopes {bi d ) of the regression of kinship coefficient values on the natural logarithm of distance {dl^ for each estuary, the coefficient of 
determination R 2 , the average and standard error (SE) of the kinship coefficient among individuals separated by less than 100 m (i.e. 0.1 km) (F A ) for each estuary, 
and the intensity of SGS (Sp) calculated for pairwise distance among individuals up to 10 Km within each estuary is shown. When GENELAND detected spatial 
genetic discontinuities within estuaries, the SGS parameters were recalculated for each genetic cluster when N > 50. Each one of the genetic clusters were named 
based on color assignation from Figure 1 including pink (P), malachite green (MG) and orange (O) within Caribbean estuaries and yellow (Y), purple (P), 
tourmaline green (TG), lapis lazuli (LL) and red (R) in Pacific estuaries. 
***P<0.001; **P<0.01; *P<0.05. 



germinans populations observed in Panama has an effect 
on the levels of outcrossing and therefore on the genetic 
structure in this species. In isolated individuals or low 
density populations geitonogamy could be an advanta- 
geous breeding system [36]. Moreover, the r=-0.64 esti- 
mates of pollen vs. seed flow suggests that although seed 
dispersal is likely equivalent to pollen dispersal in A. ger- 
minans, in a patchy matrix, this seed dispersal still is 
usually very local, leading to greater biparental inbreed- 
ing, and increasing population genetic differentiation 
and spatial genetic structure in nuclear genomes [38,39]. 

Our data also suggest that historically ambophilous 
pollen dispersal mechanism of R. mangle has been more 
efficient promoting outcrossing and long-distance gene 
flow than entomophilous pollination system of A. germi- 
nans. Based on previous genetic studies, it is predicted 
to have higher dispersal potential and therefore low gen- 
etic structure in species with wind pollination system 
over species with insect pollinator system [33]. However, 



reproductive biology studies in other ambophilous spe- 
cies suggest a completely opposite trend indicating that 
in self-compatible species as is the case of R. mangle, 
wind is actually the mechanism that promotes selfing 
and that out-crossing is associated with insect pollen 
distribution. The reason is that abiotic mechisms as 
wind do not target distant receptive flowers as efficiently 
as insects [40,41]. 

Although we do not know the rates of wind-to-insect 
pollen dispersal in R. mangle, the ratios of pollen-to-seed 
dispersal of r=7.7 estimated in this study are lower than 
those estimated in exclusively wind-dispersed plants (r=17 
[42] and r=200 [43]). In addition, based on the Hamilton 
and Millers method calculations, pollen vs. seed was not 
significantly different from each other meaning that 
ambophilous pollen dispersal is less efficient than exclu- 
sively anemophilous pollen dispersal and/or that seed dis- 
persal in this species is comparatively higher than in other 
exclusively wind pollinated species. 
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The contrast of seed viability between the two species 
may also have a strong influence in levels of genetic 
structure. Rhizophora mangle species have the highest 
longevity seeds of any mangrove genus [15,29]. Direct 
experiments regarding establishment of seeds after long 
periods of floating exposure in sea water have showed a 
60% success rate after 247 days floating for R. mangle, 
higher than any other Rhizophoracea species [44]. There 
is no similar quantitative data on establishment success 
after dispersal for A. germinans; however, its seed lon- 
gevity is shorter than R. mangle [29]. The local genetic 
structure observed in Panama largely corroborates this 
pattern, especially in the Caribbean, where two estuaries 
separated by 300 km showed identical chloroplast and 
nuclear genetic diversity in R. mangle but strong struc- 
ture in A. germinans. Although A. germinans showed 
evidence of shared cpDNA haplotypes among estuaries 
on the same coast, there are also some cpDNA haplo- 
types and nuclear alleles that are restricted to each estu- 
ary, generating strong structure, even within estuaries. 
Thus, in a variable estuarine environment where seed 
movement is stochastic [18], our data suggests that 
higher propagule longevity leads to a greater chance of 
successful establishing at long distances, increasing gene 
flow and decreasing population structure [45,46]. 

The historical imprint on genetic structure 

The Isthmus of Panama represents a 20 My to three My 
old barrier to seed gene flow between the Atlantic and 
Pacific oceans, and is the narrowest terrestrial area in 
the New World that separates mangrove populations in 
each ocean [47,48]. Based on microsatellite and chloro- 
plast genetic diversity observed across the Isthmus, this 
land mass has created high levels of genetic structure 
and a strong barrier for contemporary seed gene flow. 
Moreover, our results indicate that the Isthmus also 
represents a strong barrier to pollen flow in each species. 
This is likely due to the absence of a continuous terres- 
trial population that spans the entire distance between 
coasts. Thus, the strong isolating effect that the rise of 
the CAI has had on population differentiation for these 
two and other species highlights the role of restricted 
seed dispersal in creating spatial genetic structure in 
hydrochorious marine species [49-51]. 

Population differentiation observed between different 
oceans is exceptionally high compared to other tropical 
tree species. For example, Dick and Huertz [52] report 
an average F ST = 0.14 for microsatellites variation 
from the Neotropical tree Symphonia globulifera. Within 
Panama, S. globulifera averaged F S t = 0.11 for samples 
taken across the Isthmus of Panama, which is within the 
range of observations that we see when comparisons are 
made within oceans for the mangrove species here. How- 
ever, even trans-Andean F ST between Mesoamerican 



populations of bird pollinated, mammal dispersal S. globuli- 
fera and those in the Amazon separated by > 3,000 km 
showed maximum F ST = 0.27, still lower than lowest pair- 
wise F S t for A. germinans made between trans-Isthmian 
populations of Montijo Gulf and Costa Arriba (F S t = 0.31) 
separated by < 100 Ion. Furthermore, populations of insect pol- 
linated and wind dispersed mahogany (Sweitenia macrophyla) 
across Mesoamerica also show a lower overall F S t = 0.10, 
with the largest pairwise differences (F ST = 0.238) observed 
between Panamanian and Guatemalan populations at a dis- 
tance of > 1,600 Ion [53]. 

In spite of similar effects of CAI in the genetic struc- 
ture of these two mangroves, we found unexpectedly 
high levels of cpDNA diversity and structure in A. germi- 
nans that suggests a level of diversity that could be more 
influenced by population history and demography than 
current gene flow [54]. The cpDNA diversity observed 
in A. germinans is remarkable given the small sample 
size and short geographic distances separating the popu- 
lations both within and between coasts. This high level 
of diversity could be indicative of historical processes de- 
termining spatial genetic structure of populations. Al- 
though the fossil record of mangroves in Pleistocene is 
scarce and therefore the reconstruction of current man- 
grove distribution is very speculative, several lines of evi- 
dence suggest that mangrove ecosystem were under 
episodic crises during the Quaternary, specially asso- 
ciated to sea-level and temperature/humidity fluctua- 
tions [55-57]. In particular, it is possible that current 
populations of A. germinans represent remnants of refu- 
gial populations created during the Pleistocene [58]. In 
fact, Panamanian populations are genetically diverse 
compared to other regions, and it has been suggested that 
both ancient introgressive hybridization and secondary 
contact between A. germinans and its sister species A. 
bicolor has occurred in Panama (especially on the Paci- 
fic side), generating a hotspot of genetic diversity [31]. 

Our results in A. germinans contrast strongly with R. 
mangle, where very little cpDNA diversity was observed. 
These two species differ in their density and distribution, 
with R. mangle forming extremely dense, continuous for- 
ests near to shore, and A. germinans forming patchily 
distributed, lower-density stands at low and middle 
intertidal zones. One hypothesis that could explain the 
current distribution of cpDNA diversity is that mangrove 
populations represent relicts of much larger ancestral 
populations but that after Pleistocene-Holocene sea-level 
fluctuations, the mangrove composition shifted to a R. 
mangle -dominated community [55,57,59]. Under this 
scenario, R. mangle could have resulted in a more efficient 
colonizer than A. germinans follow a stepping-stone disper- 
sion pattern. This process combined with self-fertilization 
observed in R. mangle could be the responsible of the 
current continous and dense populations and the low 
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cpDNA and nuclear diversity observed in this species [57]. 
Alternatively, it is also possible that only certain older 
lineages of R. mangle that are well adapted to current con- 
ditions survived in the Pleistocene-Holocene sea-level fluc- 
tiations and that only those exclusive lineages recolonized 
available habitats during the Holocene or that longer time 
of R. mangle presence in neotropics generated more lost of 
diversity via genetic drift than younger A. germinans. 

Regardless of the explanation, A. germinans joins the 
ranks of tropical tree species in Panama that show com- 
plex biogeographic history with disproportionately high 
levels of cpDNA diversity relative to other parts of their 
range [60]. The historical complexity of Panama was also 
evident when we analyzed the geographical variation in 
structure with greater population differentiation in Paci- 
fic versus Caribbean estuaries for both species in nuclear 
and plastid genetic markers. In the case of R. mangle, 
introgressive hybridization between R. mangle and its 
sister species, R. racemosa, is a novel source of genetic 
variation exclusive to the Pacific [13]. Thus, the levels of 
genetic structures and patterns of biodiversity in R. 
mangle are completely different between Caribbean and 
Pacific due to independent historical processes in the oc- 
currence and sympatry of R. mangle and R. racemosa. 

The influence of coastal morphology 

Hydrochory per se is assumed to be one of the most effi- 
cient mechanisms for LDD in plants. Therefore, it is 
expected that hydrochorious plant species should have 
high gene flow among populations and low genetic 
structure, especially in local geographic areas [61-63]. 
Our data contradict this hypothesis because the two 
mangrove species resulted genetically structured among 
and within estuaries, indicating local restrictions to both 
pollen and seed dispersal, especially in A. germinans spe- 
cies. Currently, one of the major threats in mangroves is 
fragmentation and sea-level changes associated with cli- 
mate changes [10,11]. Historically, both mangrove spe- 
cies have experimented sea-level fluctuations at several 
times and climate changes [3,4,55,56], thus, our results 
suggest that historically R. mangle have maintained more 
gene dispersal of both pollen and seed dispersal than A. 
germinans. However, our results also suggested that geo- 
graphic location is important in predicting levels of gen- 
etic structure in mangroves. Pacific populations proved 
to be more structured than Caribbean populations in 
both mangrove species. One possible explanation is bio- 
logical. The Pacific coast has been characterized by an- 
cient hybridization between A. germinans and A. bicolor, 
but there is also the current scenario of introgressive 
hybridization between R. mangle and R. racemosa. Both 
hybridization processes are apparently complicating current 
levels of structure compared with Caribbean populations 
[13,31]. The other possible explanation is abiotic, including 



basin geomorphology and connectivity due to ocean cur- 
rents. Our results suggest that current or historical land- 
scape characteristics of Pacific estuaries are in some way 
enhancing pollen and seed dispersal limitations, generating 
more structure compared to Caribbean estuaries. Similarly, 
the density of A. germinans is very variable and in some 
places, for example in French Guiana A. germinans could 
be more dense and extended than R. mangle [64]. In conse- 
quence, patterns of genetic structure among and within es- 
tuaries in that region could be completely different to 
observed in Panama. Thus, although life-history traits are 
important to predict expected genetic structure, landscape 
settings are generating a variety of situations on local scales 
that complicate any prediction in terms of expected levels 
of genetic structure [20]. Therefore, although long distance 
gene flow between South America and West Africa has 
been observed in both A. germinans and R. mangle [12,13], 
estuarine geomorphology and ocean currents in Panama, 
especially on the Pacific side, seem to be more complex, 
preventing pollen and seed dispersal. 

Methods 

Study sites and sampling strategy 

We collected leaf samples from A. germinans and 
R. mangle trees within the four largest estuaries in 
Panama: Bocas del Toro and Costa Arriba in the 
Caribbean and Montijo Gulf and San Miguel Gulf in 
the Pacific (Figure 1). Within each estuary, we selected 
ten equidistant sites across the geographic contour for 
sampling. The sampling sites included rivers, streams, 
channels and shoreline areas, each with riverine or 
fringe mangrove forests. Within each sampling site we 
established transects parallel to the water's edge, and 
randomly selected a maximum of 30 adult trees for 
each mangrove species with a dbh > 10 cm and spaced 
with a minimum distance of 5 m between sampled 
trees. All selected trees were geo-referenced using a 
GPS or a compass and measuring tape. 

Microsatellite analysis 

We used 11 and 6 microsatellite loci for the analysis of 980 
A. germinans individuals and 940 R. mangle individuals re- 
spectively, following established protocols [13,65-67] but 
with modifications [68] (See Additional file 4 for details). 
We calculated the average number of alleles per locus, 
observed heterozygosity (H 0 ), expected heterozygosity (H e ), 
and the fixation index (F IS ) for each species across the four 
estuaries using GENALEX 6.0 [69]. In addition, we calcu- 
lated the outcrossing rate by hand as 1= F IS I1 + F IS [70]. 
For each locus, we tested deviations from Hardy- Weinberg 
equilibrium (HWE) and linkage equilibrium (LE) with 
GENEPOP 3.4 [71]. Also, we tested for the presence of null 
alleles and scoring problems associated with allelic stutter- 
ing or allelic dropout using MICROCHECKER 2.2.3 [72]. 
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This analysis found no stutter bands or genotyping 
errors across microsatellites. However, some loci showed a 
positive presence of null alleles, ranged from -0.04 to 0.26 
in A. germinans and from -0.18 to 0.29 in R. mangle 
(Additional File 4). 

We compared the hierarchical genetic structure of A 
germinans and R. mangle at three geographic levels: (1) 
Among four estuaries, (2) Between two estuaries along 
same coastline and (3) Within estuaries. For (1) among 
estuaries, we used a F ST based AMOVA [73] in ARLE- 
QUIN 3.5 [74] after 10,000 permutations. We compare 
this result with the ENA method implemented in 
FREENA [75] to correct the bias induced by the pres- 
ence of null alleles on the F S t estimation after 10,000 
replicates (Additional File 4). For levels (2) Between 
two estuaries along same coastline and (3) Within 
estuaries we used STRUCTURE 2.2 [76-78], GENE- 
LAND 3.1.4 [79,80] and, for mixed-mated R. mangle, 
INSTRUCT [81] . 

In STRUCTURE 2.2 we assumed an admixed model 
and a uniform prior probability of the number of popu- 
lations, K. All the runs were performed with 500,000 
MCMC replicates after a burn-in of 50,000 replicates. 
We used a model of correlated allele frequencies varying 
the level of structure from K = 1 to 6 populations. Ten 
independent runs were done for each value of K to gen- 
erate our estimate of the true number of demes [82]. 
Previous empirical and simulation analysis using 
STRUCTURE showed that null allele presence have a 
low effect in the accuracy of assignment tests [83,84]. 
This effect was moderate even in populations with a fre- 
quency of null alleles > 0.917 for a single locus [83]. 
Therefore we performed the STRUCTURE 2.2. analysis 
with no correction to the raw data to account for null 
alleles. Nevertheless, the GENELAND 3.1.4. analysis was 
performed using a explicit null alleles presence model as 
we explain below. 

The GENELAND 3.1.4 model assumes that population 
membership is structured across space. Thus, if this as- 
sumption is correct, the power of inferring clusters 
based upon the combination of genetic and geographical 
information increases compared with using STRUC- 
TURE alone [85,86]. However, in the case of weak 
spatial organization, the inferred structure of GENE- 
LAND is expected to be similar to STRUCTURE infer- 
ence [80]. For each run in GENELAND, we used the 
spatial D-model to calculate allele frequencies and set 
the maximum rate of the Poisson process and the max- 
imum number of nuclei (i.e. three times the total num- 
ber of individuals) according to the total number of 
individuals collected within each coast. In addition, we 
set an uncertainty attached to spatial coordinates fixed 
to 50 m. We performed ten independent runs of GENE- 
LAND allowing K to vary from 1-10 populations using 



the simultaneously uncorrelated allele frequency model, 
the spatial model, and the null allele model due to the 
presence of null alleles detected in both mangrove spe- 
cies (Additional File 4). We completed 100,000 iterations 
for each independent run, saving every 10 th iteration 
using a burn-in of 5,000 iterations. Finally, as mating 
systems may be different between the two species under 
study and selfing rates could influence deme structure in 
mixed-mating species, we used INSTRUCT to infer sim- 
ultaneously the selfing rates of R. mangle and its genetic 
structure at two levels: (2) Between two estuaries along 
same coastline and (3) Within estuaries. In INSTRUCT, 
HWE is not assumed as it is in STRUCTURE; rather the 
expected genotype frequencies are calculated based on 
selfing rates. We ran five independent chains, each chain 
having 500,000 iterations steps, 250,000 burn-in itera- 
tions, a thinning interval of 10 and assuming a Dirichlet 
process of mixture. 

Spatial genetic structure within estuaries 

We investigated spatial genetic structure (SGS) at level 
(3) within estuaries using a spatial autocorrelation ana- 
lysis [87]. For this analysis, we calculated the pairwise 
kinship coefficient between all individuals (F t j) separated 
at different distance classes following [88] up to 10 km. 
Kinship coefficients (Fij) were regressed using the loga- 
rithm of the spatial distance between individuals (Idij) 
within each one of the estuaries analyzed. Standard 
errors were assessed by jackknifing data over each locus. 
We generated uneven lags with constant number of 
individuals inside distance classes (N > 100) with > 90% 
of pairwise relationships among nearest neighbors 
included within the first interval. Using SPAGeDi 1.3 
[89], we calculated the regressions of kinship coefficients 
(Fy) vs. natural logarithm of distance classes (Idij) to pro- 
vide the regression slope (b id ). We tested the significance 
of SGS and the IBD in two dimensions using the 
observed slope (bid) of the linear regression of the kin- 
ship coefficient on the logarithm of the distance class 
against the null hypothesis H G : bi d = 0 (i.e. the overall 
absence of SGS) by comparing the observed values with 
those obtained after 10,000 permutations of individuals 
between locations. Where GENELAND detected a 
spatial deme within an estuary, this procedure was ap- 
plied to each inferred deme that was represented by at 
least 50 individuals, excluding admixed individuals. 

To compare the extent of SGS between the two man- 
grove species over the same geographic scales and across 
different estuaries, we calculated the Sp statistics on 
spatial scales of up to 10 km. The Sp statistics were cal- 
culated from the slope of the regression (b id ) of the Kin- 
ship coefficient [Fij of 88] against the logarithm of the 
distance, S p = - bi 10km /(l - F A ) where the regression 
slope, buokm i s l ess than 10 km and F A is the average 
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kinship coefficient between individuals belonging to the 
first distance class (< 100 m in all cases). This first dis- 
tance class included all pairs of neighbors [90]. In the 
cases where we found denies within estuaries we 
repeated this procedure for each one of the subpopula- 
tions if N > 50. 

Chloroplast analysis 

A subsample of 58 individuals of A. germinans and 60 
individuals of R. mangle and distributed across the four 
estuaries (i.e. 14-18 individuals per estuary) was 
sequenced using two chloroplast (cpDNA) non-coding 
regions atpl-atpH and psbJ-petA [91]. We redesigned 
primers to avoid short-repeat regions and improve se- 
quence quality in all regions except for the atpl-atpH re- 
gion in A. germinans [91]. Modified primers for the 
psbJ-petA region in A. germinans were F: AGATTGATC 
GATATCGGGTTC and R: GGAAAACCGAAACCCA 
GAG Modified primers for R. mangle analysis, PCR and 
sequencing specifications for both species were described 
previously [13]. For each mangrove species we calculated 
the haplotype diversity using DNASP 4.5 [92] . In addition 
we constructed a median joining network [93] for the 
combined cpDNA regions using NETWORK 4.5.1.0 
(fluxus-engineering.com). We calculated population struc- 
ture using ARLEQUIN 3.5 [74] at level (1) among four es- 
tuaries and (2) between two estuaries from same coastal 
line. Finally, we used SPAGEDI 1.3 [89] to investigate SGS 
at level (3) within estuaries, using a spatial autocorrelation 
analysis of cpDNA variation. For this analysis, we calcu- 
lated the pairwise kinship analogue coefficient between all 
individuals based on the genetic distances between haplo- 
types (Ni ; ) separated at different distance classes vs. nat- 
ural logarithm of distance classes (ld i; ) to provide the 
regression slope (bid) on spatial scales up to 10 km in a 
procedure similar to that used for microsatellites [89]. 

Comparison of pollen and seed migration rates 

In order to compare pollen and seeds migration rates 
across four estuaries we estimated F S t and F IS from 
microsatellite data and G S t from cpDNA data using 
SPAGEDI 1.3 [89]. These F-estimates were used to cal- 
culate the ratio (r) of pollen migration (m p ) to seed mi- 
gration (m s ) (i.e. m p /m s ) for each mangrove species 
following the equation r=[(l/F ST hipar - 1)(1 + F IS )] - 2 
(1/Gst mat - l)I(l/G S T mat ~ 1) [Eq. 5a, 43]. In addition, 
the same F-statistics based on microsatellite and cpDNA 
variations were used to test the null hypothesis m seed = 
m P oiien comparing the expected maternal F S t predicted 
from microsatellite variation (i.e. biparental F S t) vs. 
the observed maternal F ST estimated as G S t from ac- 
tual cpDNA variation [21]. The expected maternal F ST 
was calculated following F ST mat= {ci hipar F S T bi P ar)l 

\ a mat + ( a bipar ~ a mat)^ST bipaA where a mat - 2.0 and 



Ubipar change depending of outcrossing rate (t) calculated 
for each species from microsatellite data [Eq. 10, 21]. 
Under this procedure, the null hypothesis is rejected if 
confidence intervals (± 2SE) of observed and expected 
values fail to overlap [21]. 

Conclusions 

Although there is documented direct and genetic evi- 
dence for extreme LDD in mangrove species A. germi- 
nans and R. mangle, our data across estuaries in Panama 
showed restricted gene dispersal overall and equivalent 
rates of seed and pollen dispersal in both mangrove spe- 
cies. Rh izophora mangle showed lower gene diversity 
and lower genetic structure than A. germinans. This sug- 
gest that an amphophilous pollen syndrome combined 
with a higher propagule longevity leads to a greater 
chance of successful establishing at long distances, in- 
creasing gene flow and decreasing gene diversity and 
population structure. However, species density, coastal 
geomorphology as well as ocean currents could vary 
across ocean basins and estuaries, generating a variety of 
situations on local scales that complicate any prediction 
in terms of expected levels of genetic structure. 

Additional files 



Additional file 1: Genetic diversity of nuclear (i.e. microsatellites) 
and chloroplast genomes for Avicennia germinans (black mangrove) 
and Rhizophora mangle (red mangrove) in Panama. The number of 
individuals analyzed per estuary, the number of alleles, observed 
heterocigosity (Ho), and the expected heterocigosity (He) of 
microsatellites calculated with GENALEX 6.0 is shown. In addition, the 
number of individuals analyzed per estuary, the number of haplotypes 
and the haplotype diversity found in chloroplast using DNASP 4.5 and 
NETWORK 4.5.1.0 (fluxus-engineering.com), and the GenBank accession 
number of each individual analyzed is indicated. 

Additional file 2: Median joining network indicating and 
geographic distribution of cpDNA haplotypes found in Avicennia 
germinans (Black mangrove). Within the network, the haplotype name 
and the number of individuals per each haplotype is indicated. In 
addition, for each estuary, the geographic distribution of haplotypes and 
their frequency (i.e. pie) is indicated. 

Additional file 3: Median joining network indicating and 
geographic distribution of cpDNA haplotypes found in Rhizophora 
mangle (Red mangrove). Within the network the haplotype name and 
the number of individuals per each haplotype is indicated. In addition, 
for each estuary the geographic distribution of haplotypes and their 
frequency (i.e. pie) is indicated. 

Additional file 4: PCR conditions for 1 1 microsatellite loci for A. 
germinans (black mangrove) [65,66] and six microsatellite loci for R. 
mangle (red mangrove) [67] following established protocols [1 3,65,66] 
with modifications including three primers in the PCR thus: a dye 
tagged M13 universal forward primer (5'-CACGACGTTGTAAAACGAC-3') 
[68], primer Forward (F) and primer Reverse (R) where either F or R 
primer (indicated with asterisk *) had a tail at the 5' end that was 
identical to the Ml 3 universal forward primer sequence. Amplified 
fragments from both mangrove species were electrophoretically separated on 
ABI 3130x/ and analyzed using ABI PRISM® GeneMapper™ software version 
3.7. For each estuary is indicated the frequency of null alleles calculated by the 
Brookfield method 1 in MICROCHECKER 2.2.3 [72], 
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Additional file 5: Bayesian genetic assignment of Avicennia 
germinans (black mangrove) from two Caribbean (Bocas del Toro 
and Costa Arriba) and two Pacific (Montijo Gulf and San Miguel 
Gulf) estuaries in Panama based on STRUCTURE ver. 2.2 and 
GENELAND ver. 2.0.12. The true K for each procedure after simulations 
is indicated. 

Additional file 6: Bayesian genetic assignment of Rhizophora 
mangle (Red mangrove) from two Caribbean (Bocas del Toro and 
Costa Arriba) and two Pacific (Montijo Gulf and San Miguel Gulf) 
estuaries in Panama based on STRUCTURE ver. 2.2, INSTRUCT and 
GENELAND. 2.0.12. The true K for each procedure after simulations is 
indicated. In addition, INSTRUCT was used to help to simultaneously infer 
the selfing rates of this mixed mating species and the demic structure on 
both sides of the Isthmus. 
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